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Abstract The status and accuracy of the precision Monte Carlo generators used for luminosity measurements 
at flavour factories is reviewed. It is shown that, thanks to a considerable, long-term effort in tuned comparisons 
between the predictions of independent programs, as well as in the validation of the generators against the 
presently available calculations of the next-to-next-to-leading order QED corrections to Bhabha scattering, the 
theoretical accuracy reached by the most precise tools is of about one per mille. This error estimate is valid 
for realistic experimental cuts, appears to be quite robust and is already sufficient for very accurate luminosity 
measurements. However, recent progress and possible advances to further improve it are also discussed. 
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1 Introduction 

Precision tests of the Standard Model, indirect 
. bounds on the Higgs boson mass and constraints on 
candidate models of New Physics rely on the predic- 
tive power of quantum corrections to accurately mea- 
sured observables, such as the electroweak parameters 
and the anomalous magnetic moment of the muon [l[ . 
However, this predictability is somewhat obscured by 
the non-perturbative strong interaction effects affect- 
ing the calculation of the light quark loop contribu- 
• tion to the muon anomaly and the running of the 
electromagnetic coupling constant Q. As it is well 
known, the most satisfactory solution to this problem 
consists in rephrasing, via dispersion relations, these 
non-perturbative corrections in terms of the experi- 
mental cross section of electron-positron (e'^e") anni- 
hilation into hadrons measured in the region of low- 
lying hadronic resonances. This measurement, with 
an accuracy of the order of one per cent or even bet- 
ter, is an important goal of the physics program of the 
high-luminosity e'^e~ colliders of intermediately high 



energy, i.e. the flavour factories VEPP-2M/VEPP- 
2000, DA$NE, BEPC, CESR, KERB and PEP-II 
On the other hand, the measurement of the hadron 
production cross section, either through direct scan 
or the method of radiative return, requires, quite gen- 
erally, a detailed knowledge of the collider luminos- 
ity Q At e+e~ accelerators, the most precise and 
widely used strategy to determine the luminosity is 
by means of the relation JCdt — N/ea, where N is 
the number of events of some chosen reference pro- 
cess, e the experimental selection efficiency and a the 
theoretical cross section of the normalization process. 
Therefore, the luminosity monitoring processes must 
be characterized by high statistics, clean topology and 
small amount of background, and should be calcula- 
ble with high theoretical accuracy. This is the reason 
why the QED processes of Bhabha scattering, two 
photon and muon pair production are the normaliza- 
tion processes used at flavour factories @-@| , Bhabha 
scattering being of primary importance because of its 
large cross section. At all the flavour factories, the 
final-state particles of the above processes are selected 
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*For completeness, it is worth mentioning that a clever alternative to the standard luminosity-based measurement of the 
hadronic cross section is adopted by the BaBar collaboration [J] and is under way by the KLOE experiment 0,0]. It consists in 
normalizing the hadron production cross section measured through radiative return with radiative muon events in each bin, thus 
avoiding or, better, minimizing theoretical input and allowing for the cancellation of many common systematic uncertainties, at a 
price, however, of an increased statistical error. 
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at large scattering angles, using the same detectors 
employed for the measurement of the hadronic cross 
section. Typical experimental luminosity errors lie in 
the range between few per mille and one per cent, 
e.g. 0.3%, 0.7% and ~ 1% for present Bhabha mea- 
surements at KLOE, BaBar and CLEO, respectively, 
and at the level of ^ 0.2 0.3 for the forthcoming 
CMD-3 measurement at the VEPP-2000 0. These 
experimental requirements necessitate theoretical cal- 
culations of Bhabha scattering and, more generally, of 
normalization processes with a precision of about one 
per mille, not to spoil the total accuracy in the lumi- 
nosity measurement given by the sum in quadrature 
of the fractional experimental and theoretical uncer- 
tainty. Thanks to the effort of various groups, this 
high theoretical precision has been now successfully 
achieved, as discussed in the following. 

The main focus of the present contribution is on 
large-angle Bhabha scattering, which is the process of 
primary and wider interest and for which the largest 
part of the most up-to-date results is available. A 
complete inventory of the generators not addressed 
here and used at flavour factories for the simulation 
of QED and other processes can be found in Q . 

2 Status and theory of the most pre- 
cise luminosity generators 



(NLO) QED corrections supplemented with leading 
logarithmic contributions related to multiple soft and 
coUinear photon emission. Such ingredients, together 
with the vacuum polarization correction, are strictly 
necessary to achieve a theoretical precision down to 
the per mille level, as extensively proved in the lit- 
erature Indeed, when considering Bhabha 
scattering with typical selection cuts and at centre 
of mass (cm.) energies from (j) to B factories, the 
NLO photonic corrections amount to about 15-^20% 
(with non-logarithmically enhanced contributions at 
the few per mille level), vacuum polarization intro- 
duces a correction of several per cent and higher-order 
leading effects lie between 1-^2%. The common the- 
oretical feature of BabaYaga@NLO, BHWIDE and 
MCGPJ is the appropriate combination (matching) of 
traditional techniques for the calculation of NLO cor- 
rections (Feynman diagrams) with universal methods 
for the treatment of higher-order contributions (QED 
coUinear Structure Functions in MCGPJ, QED Par- 
ton Shower in BabaYaga@NLO and YES exponenti- 
ation in BHWIDE). Consequently, the main source 
of theoretical uncertainty of these precision genera- 
tors is due to the incomplete or approximate inclusion 
of NNLO QED corrections, whose calculation comes 
therefore to play the role of benchmark for the accu- 
racy of the programs. 



During the last few years a remarkable progress 
occurred in reducing the theoretical error contribu- 
tion to the luminosity measurements through the re- 
alization of newly conceived Monte Carlo (MC) pro- 
grams or the validation of codes used in the 90s by 
LEP/SLC experiments from the Z peak region down 
to the energy range of interest for flavour factories. 
In parallel, an impressive effort was put in the calcu- 
lation of the next-to-next-to-leading order (NNLO) 
QED corrections to Bhabha scattering (see [3| for a 
recent review), that turned out to be (and will pre- 
sumably continue to be) essential to assess the accu- 
racy of the programs. 

Two dedicated eventgenerators, BabaYaga@NLO 
l3-Eil and MCGPJ El Q, were released in 



2005-2006 to provide predictions for the cross sec- 
tion of the large-angle Bhabha process, as well as of 
two photon and muon pair production, with a the- 
oretical accuracy at the level of 0.1%. In addition, 
codes well-known since the time of LEP/SLC oper- 
ation, like BHWIDE 17|, were extensively used by 



the experimentalists in data analysis at flavour fac- 
tories. All these programs include, albeit according 
to different formulations, exact next-to-leading order 



3 Accuracy of the luminosity genera- 
tors vs. NNLO QED calculations 



It is customary to define the total accuracy of the 
generators in terms of a so-called technical and the- 
oretical (or physical) precision. As independent pro- 
grams implement, as discussed above, radiative cor- 
rections of the same physical origin but according to 
different theoretical and computational details, this 
may give rise to possible differences between their 
predictions. This defines what is known as techni- 
cal precision. However, the most important source of 
uncertainty is that induced by missing or partially ac- 
counted contributions, like e.g. NNLO terms present 
in the exact perturbative calculations. These miss- 
ing or approximately included ingredients define the 
so-called theoretical accuracy. During the Working 
Group (WG) on "Radiative corrections and MC gen- 
erators for low energies" (sj a significant work was 
done to deeply investigate the size of both uncertain- 
ties. The results summarized in the following are the 
main conclusions of this effort. 
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3.1 Technical precision 

The typical strategy to settle the technical pre- 
cision of theoretical tools was initiated at CERN in 
the 90s during the WGs on precision physics at LEP 
and consists in performing tuned comparisons be- 
tween the predictions of the codes independently de- 
veloped by different groups. The word tuned means 
that the comparisons, to be meaningful, must be per- 
formed using the same set of parameters and cuts, 
in order to really test the reliability of the programs 
in all the technical details underlying their rather 
complex structure. Several tests of this type for the 
Bhabha process, according to different selection cuts 
and cm. energies, are documented in Q. It emerges 
that BabaYaga@NLO, BHWIDE and MCGPJ typ- 
ically agree within ~ 0.1% for integrated cross sec- 
tions and some per mille for distributions. There- 
fore, the work on tuned comparisons showed that the 
technical precision of the MC luminosity generators 
is well under control, the remaining differences being 
due to understood different details in the implemen- 
tation of theoretical ingredients featuring the same 
physics (factorized vs. additive formulations, differ- 
ent realization of exponentiation effects). 

3.2 Theoretical accuracy 

Concerning the theoretical accuracy of the gen- 
erators, there are essentially two possible procedures 
to assess it. The first one is to compare the pre- 
dictions of exact NNLO QED calculations with the 
corresponding O(a^) results of the generators, if the 
above perturbative calculations are available in the 
literature, as it is presently true for Bhabha scatter- 
ing. The second one is to estimate the size of partially 
accounted higher-order corrections, if exact or com- 
plete calculations are not yet available, or, as it would 
better to say now, were not available till recently. 
Actually, as further discussed in the following, the 
remaining ingredients necessary for complete NNLO 
predictions, unavailable until a few months ago, are 
now at hand: i) the full calculation of lepton and pion 
pair corrections, accounting for real pair emission be- 
yond the soft approximation [18[; ii) the calculation of 
the complete set of one-loop corrections to the radia- 
tive process e+e~ ^7*(s),7*(t) — >■ e+e~7 [19|]. How- 
ever, these ultimate ingredients appeared in the lit- 
erature just a few weeks before the beginning of the 
conference, a too limited time to perform all the phe- 
nomenological investigations necessary to obtain and 
present an update of the error estimate based also on 
the inclusion of these new calculations. 



Considering first those NNLO corrections already 
under control thanks to independent calculations and 
rather deeply investigated in connection with the pre- 
cision of the generators, they can be divided in three 
gauge-invariant classes: 

1. photonic corrections, computed in [2(?| in the 
soft-photon approximation; 

2. Electron loop corrections, calculated indepen- 
dently by two groups in 2l|, |22| and found to 
be in perfect agreement; 

3. Heavy fermion and hadron loop corrections due 
to the insertion of /i, r leptons and hadrons in 



the photon propagators, computed in 2^2 

As shown by these calculations and summarized in Q , 
the typical size of the above NNLO corrections to the 
Bhabha differential cross section is of the order of 
some per mille as a function of the electron scatter- 
ing angle, both at 4> and B factories. Not surpris- 
ingly, they are largely dominated by photonic cor- 
rections and, to a less extent, by electron loop con- 
tributions. However, the bulk of such corrections is 
included in the generators as a bonus of the matching 
procedure, that allows to effectively incorporate the 
most important piece of the NNLO photonic correc- 
tions, given by the infrared-enhanced 0{a^L) terms 
[irl {L being the collinear logarihtm). Moreover, also 
the dominant part of the electron, heavy fermion and 
hadron NNLO corrections, originating from the re- 
ducible two- loop diagrams fiil, is included at O(a^), 
albeit partially, through the inclusion of the vacuum 
polarization effects in the NLO contribution of the 
matched formulation of the generators. Hence, the 
relevant question is to what extent the exact NNLO 
results compare with the predictions of the codes at 
the same perturbative order. 

To answer this question, detailed comparisons 
were carried out between the calculations 20|, |21| 



and the corresponding results of BabaYaga@NLO. It 
turns out that, as a consequence of the factorization 
property of the infrared and collinear singularities, as 
well as of the dominance of two-loop reducible graphs 
within the full set of NNLO contributions, the relative 
deviation between the BabaYaga@NLO predictions 
and the results of 2^, 21| does not exceed 2x10"'', 
an order of magnitude below the target one per mille 
precision. An even smaller, completely negligible dis- 
crepancy, at the level of lxlO~^, was registered be- 
tween the BabaYaga@NLO calculation of the two 
hard photon emission process e"'"e~ — e+e~77 (a fur- 
ther ingredient of the complete NNLO Bhabha cross 
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section) and the exact diagrammatic answer [l3| . 

Concerning those NNLO corrections only prelimi- 
narly scrutinized and requiring further work with re- 
lation to the accuracy of the generators, they are, as 
already remarked: 

1. the exact NLO soft plus virtual QED correc- 
tions to hard photon emission in full (s- plus 
i-channel) Bhabha scattering; 

2. the complete lepton and hadron pair correc- 
tions, containing, in addition to the exact vir- 
tual and real soft-photon contributions, also 
the exact four-particle matrix elements, like 
e+e" e^e~ = e,^,T, and e+e" 
e^e~ {tt^tt~). The parentheses mean that those 
two final-state particles must be undetected to 
contribute to the Bhabha signature. 

The exact calculation of the one-loop corrections 
to e+e" — e+e~7 was performed recently in [l9l |. 
also thanks to innovative advances in the reduction 
of complex loop structures. A corresponding numer- 
ical code to derive results of phenomenological inter- 
est (for /i+/x^7 production too) is available as well. 
However, in the absence as of today of explicit com- 
parisons between the approximate results of the gen- 
erators (again, possible as a by-product of the match- 
ing) and the exact calculation, the present estimate 
of the theoretical uncertainty in the calculation of 
this NNLO ingredient must necessarily rely on par- 
tial results known in the literature 28-3^. The lat- 
ter analyses apply to t- and s-channel processes sep- 
arately and neglect contributions, like pentagon-like 
diagrams, taken into account in the exact one-loop 
calculation [l^. However, from the studies (isl - liol 
one can argue that the error induced by the NLO 
corrections to the radiative Bhabha process conser- 
vatively amounts to about 5 x 10^*. The uncer- 
tainty coming from the contribution of lepton and 
hadron pairs completes the picture of the sources of 
theoretical error. In this context, a DESY Zeuthen 
and Katowice theory collaboration [3] recently per- 
formed a complete calculation of such contribution 
using all the necessary virtual and real pair ingre- 
dients discussed above. The approximate results of 
BabaYaga@NLO, that include part of the virtual pair 
corrections thanks to the insertion of the vacuum po- 
larization in the NLO diagrams but does not take into 
account real pair effects and the two-loop pair contri- 
butions to the electron form factor, were preliminary 
compared in [s^ (for lepton pairs only) with the exact 
DESY Zeuthen-Katowice calculation. It comes out 
that the approximation of BabaYaga@NLO agrees 



rather well with the complete prediction, showing rel- 
ative deviations not exceeding few units in 10^* in the 
presence of realistic luminosity cuts for KLOE and 
BaBar experiments. 

A final, parametric contribution to the theory 
uncertainty comes from the non-perturbative light 
quark contribution Aal^^^^{q'^) to the vacuum polar- 
ization. It can be estimated by using the routines, like 
HADR5N [ail and HMNT that parametrize 

this correction in terms of the low-energy hadronic 
cross section. Indeed these parameterizations return, 
in addition to Aa'^^^^{q^), the experimentally-driven 
error Jhadr on its value. Therefore an estimate of this 
uncertainty can be simply obtained by computing the 
Bhabha cross section with Aa|j^'j^(g^)±5hadr and tak- 
ing the difference of the predictions as a measure of 
the associated theoretical error. One concludes that 
the spread between the minimum/ maximum values 
and the central ones as returned by the two routines 
well agrees and does not exceed 2xl0~* at (p factories 
and 3x10"'' at B factories [sf. This can be simply 
understood in terms of the dominance of t-channel 
photon exchange for large-angle Bhabha scattering at 
fiavour factories and the excellent agreement between 
the two parameterizations for space-like virtualities. 

The total theoretical uncertainty of the most pre- 
cise luminosity codes can be obtained by summing 
linearly the different sources of error discussed above, 
as shown in Table [TJ As can be seen, the total the- 
oretical error amounts to about 0.1% at (/), r-charm 
and B factories ^3,]. 

Table 1. Summary of the different sources of 
theoretical uncertainty for the most precise 
generators used for luminosity measurements 
at flavour factories and the corresponding to- 
tal theoretical error for the calculation of the 
large-angle Bhabha cross section. From f3|. 



Source of error (%) 





T-charm 


B 




0.00 


0.01 


0.03 


[32] 


0.02 


0.01 


0.02 


l^ivl 


0.02 


0.02 


0.02 




0.00 


0.00 


0.00 




0.05 


0.05 


0.05 


|(5<"^": 1 

1 pairs 1 


0.05 


0.1 


0.02 


ircrr | 
1 "total 1 


0.12H-0.14 


0.18 


O.llH-0.12 



4 Conclusions 

Stimulated by important experimental progress in 
the measurement of the low-energy hadronic cross 
section worldwide, the theoretical contribution to the 
uncertainty in the luminosity monitoring at fiavour 
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factories was significantly reduced in tlie last years 
down to the one per mille level. 

From one side, the development of new gen- 
erators for the simulation of the Bhabha process 
(BabaYaga@NLO and MCGPJ) and the adoption of 
the LEP/SLC tool BHWIDE paved the way to theo- 
retically consistent comparisons between the predic- 
tions of these codes for all the observables measured 
by the experiments. The consistency stems from the 
fact that all these three generators include, albeit ac- 
cording to different details, the strictly necessary the- 
oretical ingredients given by NLO QED corrections 
matched with resummation, together with the con- 
tribution of vacuum polarization. Extensive work on 
tuned comparisons between the results of the lumi- 
nosity tools greatly benefited of the activity of the 
WG on "Radiative corrections and MC generators for 
low energies" [1] and led to the conclusion that the 
large-angle Bhabha predictions of BabaYaga@NLO, 
BHWIDE and MCGPJ typically agree within ~ 0.1% 
for integrated cross sections and some per mille for 
distributions. These results, that take into account 
realistic selection criteria, demonstrate that the tech- 
nical precision of these programs in their implemen- 
tation of the radiative corrections is accurately under 
control, with a precision definitely sufficient in com- 
parison with the present experimental errors. Actu- 
ally, all these three codes are employed by the ex- 
perimentalists to precisely monitor the luminosity of 
flavour factories, generally using more than one tool. 

From the other side, various detailed cross-checks 
between the results of the generators and those of the 
presently available calculations of the NNLO QED 
corrections to the Bhabha cross section allowed to 
assess the theoretical accuracy of the MC tools on 
solid grounds. Indeed, the luminosity generators in- 
clude, by construction, contributions of 0{a^) only 
approximately. Therefore, the comparison between 
the O(a^) predictions of the generators and those 
of the exact perturbative calculations provides a re- 
liable evaluation of the theoretical accuracy of the 
codes. The sum of the various sources of uncer- 
tainty (e.g. NNLO photonic contributions, lepton 
and hadron pair corrections, one-loop corrections to 
hard bremsstrahlung) is presently at the level of one 
per mille, again when considering selection cuts of 
actual experimental interest. This error estimate ap- 
pears to be quite robust and is presently dominated, 
as shown in Table [U by the uncertainties associated 
to the MC calculation of pair corrections (amount- 



ing to a few units in 10"'') and of e+e~ — e+e~7 
with NLO accuracy (at the level of ^ 5 x 10"''). It 
is comparable to that achieved about a decade ago 
in the luminosity measurement at LEP/SLC collid- 
ers through small-angle Bhabha scattering [33]. 

In spite of this remarkable progress, further ad- 
vances would be desirable. For example, the work on 
tuned comparisons was mainly made for the Bhabha 
process. It would be worthwhile to put forward such 
studies to cover the other luminosity reactions of in- 
terest, by comparing, for instance, the precision pre- 
dictions of BabaYaga@NLO and MCGPJ for two pho- 
ton production 0, Q and the results of the whole set 
of the available programs for the simulation of the 
processes of /i'*'/^" and /i+//~7 production. For the 
latter processes some results can be found in 0, . 
A better assessment of the theoretical accuracy could 
be achieved along the following directions. The pre- 
liminary analysis of the effect of lepton pair correc- 
tions performed in ^ should be extended to the study 
of hadron pairs and to a more thorough investigation 
of the pair correction dependence from the cuts of ex- 
perimental interest. On this topic, work is already in 
a rather good shape but still in progress [35| . The re- 
sults of the generators for the one-loop corrections to 
the radiative Bhabha signature should be compared 
with the predictions of the now available exact calcu- 
lation [l9j , again in the presence of realistic selection 
criteria. Such an analysis would be also important 
for the normalization process e+e" //+/i"7. This 
study is ongoing but is still at a very preliminary 
stage [3^. Furthermore, following the same strategy 
used for Bhabha scattering, it would be also inter- 
esting to quantify to what extent the presently avail- 
able parameterizations of Aajj^dr('?^)i that manifest 
some discrepancies in their predictions for time-like 
momenta [2 01 , affect the uncertainty of the calcula- 
tion of the cross section of s-channel annihilation into 
/i+/i" and /i+/i~7. 

All these advances would eventually lead to im- 
provements or refinements in the luminosity codes 
and possibly allow a control of the theoretical error 
in luminosity measurements at present flavour facto- 
ries and at future super- _B colliders below the one per 
mille barrier. 

Guido Montagna thanks the organizers for the 
kind invitation and the very friendly atmosphere of 
a perfectly organized and interesting workshop. 
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